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Abstract 


A  uniform  magnetic  field  is  required  to  calibrate  magnetic  induction  sensor  coils.  Two 
types  of  coil  designs  are  considered;  a  single  finite  length  solenoid  and  a  coaxially  split 
coil  pair  (Helmholtz  coil).  The  design  must  be  able  to  reliably  predict  the  magnitude  of 
the  calibration  field  and  its  homogeneity,  both  radially  and  axially,  to  within  certain 
tolerances.  The  experimental  results  of  a  design  based  on  the  split  coil  pair  are  given. 
Furthermore,  it  is  suggested  that  the  quasi-uniform  field  in  the  mid-plane  of  a  split  coil 
pair  can  be  usefully  extended  by  narrowing  the  gap  of  the  conventional  Helmholtz  coil 
pair  spacing.  In  such  cases,  the  flux  in  a  circular  region  in  the  mid-plane  can 
approach  the  value  that  would  exist  for  a  perfectly  uniform  field  because  of  the  partial 
cancellation  between  fields  greater  than  and  less  than  the  central  field.  The  deviation 
of  the  central  field  from  that  of  an  infinitely  long  solenoid  is  also  discussed. 
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Magnetic  Field  Distribution  Within 
Uniform  Current  Density  Coils  of  Non- 
Zero  Cross-Section  and  the 
Design  of  Helmholtz  Coils 

1.  Introduction 

It  is  often  necessaiy  to  have  a  source  for  a  lutifoim  magnetic  field,  for  example, 
to  provide  a  held  of  known  uniformity  and  magnitude  for  calibration  purposes. 
A  single  coil  is  easier  to  construct  than  a  combination  of  two  or  more  coils  but 
its  field  homogeneity  may  be  inadequate;  thus  a  quantitative  comparison  is 
required.  Furthermore,  there  appears  to  be  a  paucity  of  data  and  tables 
available  in  the  literature  that  allow  the  design  of  a  system  to  generate  a 
homogeneous  field  over  a  given  region.  The  coil  geometries  represented  in  the 
available  tables  can  be  too  restrictive  and  the  given  error  limits  or  uniformity 
tolerances  may  be  unsuitable.  Another  important  feature  which  is  not 
emphasized  is  the  field  vmiformity  of  a  slightly  uncompensated  coil  pair  having 
a  gap  which  differs  from  the  optimiun  coil  pair  spacing  of  a  true  Helmholtz 
coil.  Such  a  coil  can  generate  a  nearly  imiform  field  over  a  larger  region  than 
the  true  Helmholtz  coil.  If  the  flux  of  the  axial  field  is  required,  then  the 
introduced  error  is  improved  significantly  because  of  the  possible  cancellation 
of  field  non-uniformities. 

In  this  report,  a  brief  description  of  solenoid  design  is  given  in  terms  of  the 
geometry  of  the  coU  [1].  There  are  two  reasons  for  examining  these  details. 
Firstly,  it  enables  the  axial  field  homogeneity  to  be  determined  at  radial 
distances  away  from  the  coil  centre.  Secondly,  as  pairs  of  coils  will  be 
considered  as  candidates  for  providing  homogeneous  fields,  the  formalism 
introduced  here  for  the  single  coil  will  also  apply  to  more  complex  coil  designs. 


1.1  Magnetic  Field  at  Centre  of  a  Coil  of  Non-zero 
Cross-section 

The  magnetic  field  along  the  axis  of  a  circular  current  loop  is  (Fig.  1) 


H,(2,0) 


I 


(1) 


where  H  is  in  units  of  A  a  is  the  loop  radius  in  m,  z  is  the  distance  along 
the  axis  in  m  and  I  is  the  current  (amps)  in  the  conductor.  Thus 

Ho  =  H^(0,0) 

-  J_ 

"  2a 


and 


H,(z,0) 


H„ 


(3) 


Figure  1:  Magnetic  field  components  for  a  circular  current  loop;  and  are  the 
radial  and  axial  components  respectively. 
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The  coil  geometry  is  defined  in  terms  of  two  parameters  a  and  P,  see  Figiire  2. 
In  this  report,  where  appropriate  in  mathematical  expressions,  the  following 
symbols  will  be  used  for  dimensionless  distances:  ^  =  r/a|  (6  =  n/1),  y  =  z/O] 

(0  =  0)  and  p/fl|  (0  ^  it/2,  0  ^  0).  Either  expression  (e.g.  4  or  r/a^)  will  be  used 
interchangably. 


p  =  b/a-]  ^  ^ 


Figure  2:  Definition  of  coil  geometry  parameters  a  and  ^  for  a  coil  of  non-zero  cross- 
section.  For  a  current  sheet  a^  -¥02  =  a,  and  cos  0  =  a  =  1,  p  =  b/a. 


The  current  density  is  /  (.Am'^)  (=  /^/cross-section  area)  and  for  N  turns  of 
uniform  current  density. 


M 

^  "  2b(a2-a^y 


(4) 


If  the  space  factor  X  (which  represents  the  active  cross-sectional  area  of  the 
winding  divided  by  the  total  cross-sectional  area)  is  not  unity,  then 


9 


/  - 


(5) 


where  /j  is  the  current  density  in  the  conductor.  As  A]  02,  a  current  sheet 
obtains,  i.e.  /'  =  NI/2b,  however,  now  /'  is  equivalent  to  current  per  unit  length 
and  I  oe  /'  dz.  Thus  for  a  current  shert  obtains  from  eqn.  1 


1 

NI 


(6) 


1.1.1  Non-zero  Crosssection  Coil 

For  a  non-zero  cross-section  coil  (Fig.  2),  the  held  at  the  centre  of  the  coil  Hq  is 
given  by  (see  eqn.  1) 


Hr 


=  ;^Xfl,pln 


*b 

-b 

( 


~{p~7zP^ 

a+(o^  +  P^)’^^ 


dzdr 


1  ♦(!  +p2)V2 


(7) 


The  held  Hg  is  signihcant  iri  the  following  analysis  because  the  helds  at  any 
point  along  the  axis,  or  away  from  the  axis  may  be  expressed  in  terms  of  the 
central  held  Hq.  Furthermore,  this  also  applies  to  the  case  of  multiple  concentric 
coils. 

One  may  conveniently  separate  an  entirely  geometry-dependent  held  factor 
f(a,  p)  from  eqn.  7  so  that  Hq  can  be  readily  estimated  from  values  of  f  (a,  P) 
represented  in  graphical  form  (see  Figs  3  and  4)  or  tabular  form,  and  the 
variation  in  Hq  as  a  function  of  coil  geometry  can  also  be  readily  determined. 


F(o,P)  ■  p  X  In 


>1 


(8) 


p[sinh'’(a/p)  -sinh'’(l/P)], 


(9) 


Alpha, 


Alpha, 
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Note  that 


(10) 


f (a,  -p)  -  -f (a,p). 


Now 


-^Xa,F(a,p) 

^  m  ^ 


2Pflj(a-l) 


f(a,p). 


(11) 


1.1.2  Deviation  of  a  Field  from  an  Infinitely  Long  Solenoid 

Figures  3  and  4  also  display  the  variation  in  Hg  from  the  ideal  infinitely  long 
solenoid  given  by  H  =  nl  where  n(=  N/2Pflj)  is  the  turns  per  unit  length. 

For  example,  for  a  fixed  small  value  of  a(a  =  15),  f(a,  P)  is  approximately 
constant  for  p  greater  than  2.0  and  thus  Hg  is  proportional  to  nl  for  a  coil 
having  a  dimensionless  length  of  2P  >  4.0.  Taking  this  point  further,  as  P  ->  <», 
eqn.  9  gives 


F(a,p)  =  (a-D- 


(0^-1) 

6P" 


^3(a5-l) 

40p^ 


(12) 


!t  is  clear  that  the  central  field  of  a  real  solenoid,  Eqn.  11,  approaches  that  of  an 
infinitely  long  solenoid,  when  f(a,  P)/(a  -  1)  1.0.  In  Figure  5,  the  deviation 

factor  F(a,  P)/(a  - 1)  is  plotted  as  a  function  of  the  normalized  coil  length  (2P) 
for  various  normalized  coil  thicknesses  (a  -  1).  It  is  now  easy  to  visualize  to 
what  extent  the  coil  thickness  and  length  affects  the  deviation  of  the  central 
field  from  that  of  an  infinitely  long  solenoid  (f(a,  P)/(a  - 1)  =  1.0). 


2.  Non-Central  Field 


2.1  Introduction 

The  axial  field  away  from  the  centre  of  the  coil,  n/1),  see  Figure  2,  may 

be  computed  from  published  tables  given  by  Hart  [2].  Whilst  these  are 
comprehensive,  there  are  several  disadvantages.  Interpolation  may  be  required 
(over  a  surface  or  line)  and  the  tables  do  not  easily  lend  themselves  to  the 
calculation  of  parameters  required  to  attain  a  specific  field  homogeneity,  either 
from  one  coil  or  by  a  combination  of  several  coils.  One  useful  feature  of  Hart's 
tables  is  that  the  fields  may  be  computed  at  points  well  away  from  the  central 
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region  of  the  coil,  even  at  points  within  the  cuirent  distribution,  provided  the 
cunent  flow  is  uniform. 


Figure  5:  The  deviation  j'actor  F(a,  Pl/(a  -Vasa  function  of  coil  length  2p  for 
various  coil  thicknesses  (a  -  V. 


A  superior  technique  introduced  by  Garrett  [3]  treats  the  field  of  an  axially 
symmetric  system  as  an  exf>ansion  of  zonal  harmonics  (legendre  polynomials), 
llie  significance  of  this  approach  is  that  the  held  at  any  point  in  the  vicinity  of 
the  origin  can  be  express^  as  an  infinite  series  of  terms  involving  simple 
algebraic  expressions.  The  vicinity  specifies  a  central  region  whose  boundary 
is  a  sphere  that  contains  no  field  source  and  satisfies  Laplace's  equation.  (A 
source  is  any  discontinuity  of  the  windings,  for  example,  the  boimdary  edge  of 
a  solenoid.)  Thus  the  distinctive  feature  of  the  zonal  harmonic  method  is  that 
it  permits  a  clean  separation  of  the  field  geometry  fivm  that  of  the  generating 
system  and  it  allows  the  field  to  be  described  as  a  continuous  scalar  function  of 
the  field  coordinates  [4]. 


2.2  Field  in  a  Central  Zone 

The  magnetic  field  in  the  central  zone  of  a  cylindrically  symmetric  nugnet  can 
be  written  as  a  power  series,  involving  Legendre  polynomials,  which  converges 
everywhere  within  a  sphere  excluding  any  field  source.  If  the  origin  is  on  the 
midplane  of  symmetry,  the  series  expansion  will  have  only  even  terms  and  the 
axial  and  radial  fields,  and  H,,  can  be  expressed  as  [3],  [5] 


(13) 


H,(p,0)  =  -sinOj:  <^4) 


^  (2n-l)! 


where 


H^^‘^’(0,0)  = 


(z,0)^ 


dz 


>.o 


(15) 


P„{u)  and  FJ,u)  are  respectively  the  Legendre  polynomials  and  their  first 
derivatives  with  respect  to  cos  0  («  =  cos  0),  see  Table  1,  and  p  and  0  are 
defined  in  Figure  1.  Equations  13  and  14  reduce  to  [1],  [5] 


E  E 

H,(p/a,,0)  -  HolO--i(p/fl,)^ P2(M)--i(p/fl,)^ ?*'(«) -...Isin0  <^7) 
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where  Hg  is  defined  in  eqn.  7, 


E 


2r 


1 

'd^H^{z,0) 

Ho(2n)! 

1 

(18) 


and  ^2nf  for  a  uniform  current  density  coil  [1],  are  given  in  Table  2.  Note  that 
in  eqns.  16  and  17,  the  fields  in  a  region  of  space  around  the  origin  depend  only 
on  Hq  and  a  geometry  factor. 

Table  1:  Even  order  Legendre  polynomials  and  Oieir  first  derivatwe  xuith  respect  to 
cos  d,  u  =  cos  6 


Po(u)  =  1 

P2(u)  =  1/2  (3tt2 . 1) 

P^(u)  =  1/8  (35u<  -  30«2  +  3) 

Ps(u)  =  1/16  (231«‘  -  315m*  +  lOSu^  -  5) 
Pg(u)  =  1/128  (6435m*  -  12012m‘  +  6930m* 
-  1260m^  +  35) 

M  =  1  (0  =  0) 

P,(±  1)  =  (±  D" 

(±  1)  =  (±  1)"*'  X  «(n+l)/2 


P'o(«)  =  0 

F2(u)  =  (6m)/2 

P'gtu)  =  1/8  (140m*  -  60m) 

P'tiu)  =  1/16  (1386m*  -  1260m*  +  210m) 

P'glM)  =  1/128  (51480m^  -  72072m*  +  27720m* 

-  2520m) 

M  =  0  (0  =  n/2) 

P„(0)  =  On  :  odd 

P,(0)  =  (-  1)"/*  X  [1.3.5..(n-l)/2.4.6..nl  n  ;  even 
P'g(0>  =  «P^,(0) 


As  shown  in  section  3.1,  the  axial  field  at  any  point  along  the  axis  is 


H,(y,0)  = 


fM 


(Y  +  P)ln 


-(7-P)ln 


a  +  [a2  +  (7  +  p)2r/2 


[l+ll+(Y  +  P)2p2 
a+Io^+(Y-P)^P^^ 


1*[1+(Y-P)"P"  J 


(19) 


where  y  =  z/a-^  is  the  distance  along  the  z-axis  normalized  to  the  inside  radius 
of  the  coil,  see  Figure  2.  Equation  7  obtains  when  y  =  0. 
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Table  2:  Error  confidents  for  uniform  current  density  coils 


r-P^ 


f(a,p)E,(a,P)  =  ’  [Cf'"-C3®^"l 


f(a,p)E4(a,p)  -  _Llcf^2+3C2+15C2) 
24  P 

-  (2 +30^+1504)] 


f(a,p)E,(a,p)  =  _L,(cf^8+12C,+15C2-70C2+315C2S 
240p® 

-  Ci'^(8  +I2C4  +I5C4  -7OC4  +315C4^] 


f(a,p)E8(a,p)  =  _L_[Cf^l6+24C2+30C2 +35C2+315C2 

896P^ 

-  2079C|  ♦3OO3C2) 

-  (16 *2404  +3OC4  +35C4  +315C4 

-  2079C4 +3OO3C4)] 


KdHJdz\^  =  (d%/dz^\^  =  {d^^^HJdz^\^  *  0 


2.3  Field  Uniformity  for  a  Single  Coil 

Along  the  z-axis  («  =  1,  see  Table  1),  H,  =  0  by  symmetry  and  is  given  by 
eqn.  16, 

H,(y,0)  =  Hod  +  E2(a,p)T^  +  E4(a,P)/  +  ...).  <2 


Further,  at  0  =  n/2,  the  axial  field  is 
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(21) 


-  H„[l  -^E;(a,p)5=*|E,((i,p)5‘  -^E,(o,p)5‘ 


For  a  simple  solenoid,  £2(a,  P)  is  negative  and  the  axial  field  decreases  along 
the  z-axis,  eqn.  20.  However,  the  axial  field  increases  radially  from  the  centre, 
eqn.  21.  Close  to  the  centre,  where  the  £2(a,  P)  term  predominates,  the  axial 
field  deviation  from  Hg  in  the  radial  direction,  n/2)/Ho,  is  only  half  that  of 
the  axial  direction.  Thus  eqn.  20  can  be  used  to  determine  a  bound  for  the 
homogeneity  of  the  axial  held  within  a  sphere  of  radius  y  (<  a)  since  the 
inhomogeneity  associated  with  any  individual  term  is  greatest  in  the  z  direction 
(see  also  section  35). 

It  is  now  necessary  to  consider  the  uniformity  of  the  axial  field,  in  terms  of 
coil  geometry  a  and  P,  over  a  sphere  of  given  radius  as  a  percentage  deviation 
of  the  central  field,  that  is 


H,(p/«,,e) 


*  1  +  H' 


(22) 


where  H  is  obtained  from  eqns.  20  and  21  for  0  =  0  (p/a^  =  y)  and  0  =  n/2 
(p/fli  =  %)  respectively.  The  ratio  H^fp/aj,  0)/Ho  has  been  evaluated  for  three 
values  of  p/aj;  0.30, 050  and  0.707  (1/V2).  The  results  are  sho'wn  graphically 
in  Figures  6  to  8.  The  terms  in  eqns.  20  and  21  were  computed  in  double 
precision  up  to  and  including  the  eighth  order  contribution.  This  allowed  the 
convergence  of  the  series  expansion  to  be  examined.  A  selected  number  of 
£„(a,  P)^"  (n  =  2, 4,  6,  8)  results  are  given  in  Table  3;  these  are  applicable  to 
both  H,(y,  0)/Hg  and  h/D/Hq.  It  is  worth  noting  that  the  convergence  of 
the  Hj(4,  Jt/2)/Hg  series  is  better  than  that  indicated  by  Table  3  because  the 
coefficients  1/2,  3/8, 5/16  and  35/128  were  omitted  to  allow  the  results  to 
apply  to  both  series. 

The  convergence  of  the  series  expansion  (eqns.  20  and  21)  is  such  that  for 
most  values  of  a  and  P,  only  the  first  two  terms  are  required,  depending  on  the 
desired  accuracy  although  clearly,  more  terms  must  be  included  as  (p/aj) 
approaches  unity.  However  convei^ence  is  slow  when  P  is  small,  especially 
when  %  >  05.  As  a  worst  case,  when  a  =  1.01  and  P  =  0.1,  the  eighth-order 
term  £9(0,  P)^*  contributes  respectively  0.1%  [0.05%],  3%  [0.9%]  and  21%  [5.4%] 
for  (p/fli)  =  0.3,  05  and  0.707  in  the  series  expansion  of  0)/Hq  and 

H, (^,  n/2)/Ho. 

As  shown  in  Figures  6  to  8,  for  a  given  coil  thickness  a,  the  deviation  of  the 
axial  field  from  the  central  field  decreases  as  the  coil  lengthens.  Conversely,  as 
the  coil  thickness  increases  for  a  fixed  p,  the  deviation  H^/Hg  decreases. 
Furthermore,  above  P  -  15  there  is  little  change  in  as  a  function  of  a  and 
hence  P  appears  to  play  the  dominant  role  in  determining  field  homogeneity. 
(This  case  approaches  that  of  an  infinitely  long  single  layer  solenoid,  see  section 

I. 1.2.)  For  lower  values  of  P,  the  field  homogeneity  is  sensitive  to  a.  At 
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P  =  1.0,  ^  =  0.707,  n/2)/Ho  varies  by  only  4%  fron»  1.088  at  a  =  1.01  to 

1.053  at  a  =  5.0.  For  ^  =  0.5,  the  variation  is  less  than  2%. 


0  0  0-2  OU  0  6  0-8  1  0  1  2  1-4  16  1-8  20 


Figure  6:  Deviation  of  axial  fields  HJlr/aj,  n/21  and  HJ[zlap  0)  from  the  central  fidd 
Hq  of  a  single  coil  as  a  function  of  a  and  P  for  p/aj  =  0.3. 
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Figure  8:  Deviation  of  axial  fields  H/r/aj,  n/2)  and  H/z/aj,  0)  from  the  central  field 
Hq  of  a  single  coil  as  a  function  of  a  and  P  for  p/a,  =  0.707. 


Table  3:  Selected  values  of  E„(a,  terms 


a 

P 

£2(0, 

£4(0,  p)4* 

£6(a,  P)4‘ 

Es(a,  p)4« 

1.01 

0.1 

-1.31E-1 

^  =  0.30 

1.4E-2 

-1.4E-3 

1.3E-4 

1.0 

-3.4  E-2 

-3.0E-4 

35E-5 

45E-7 

2.0 

-5.4  E-3 

-l.lE-4 

-l.lE-6 

2.7E-9 

35 

-7.8  E-4 

-7.6E-6 

-5.7E-8 

-3.4E-10 

5.0 

-2.0  E-4 

-l.lE-6 

-4.7E-9 

-1.8E-11 

2.0 

0.1 

-7.2  E-2 

4.9E-3 

-3.5E-4 

2.5E-5 

1.0 

-19  E-2 

1.6E-4 

llE-5 

-1.7E-7 

2.0 

-7.4  E-3 

-6.9E-5 

7.4E-8 

9.2E-9 

35 

-1.4  E-3 

-l.OE-5 

-4.8E-8 

-1.2E-10 

5.0 

-4.1  E-4 

-1.9E-6 

-6.4E-9 

-1.7E-11 

5.0 

0.1 

-4.0  E-2 

2.3E-3 

-1.5E-4 

l.lE-5 

1.0 

-1.9  E-2 

1.3E-4 

4.6E-6 

-6.8E-8 

2.0 

-6.9  E-3 

-1.6E-5 

1.2E-7 

2.8E-9 

35 

-12  E-3 

-5.6E-6 

-7.3E-9 

4.1E-12 

5.0 

-9.1  E-4 

-1.9E-6 

-14E-9 

-1.9E-12 

1.01 

0.1 

.3.6  E-1 

^  =  0.50 

l.lE-1 

-3.0E-2 

7.9E-3 

1.0 

-9.4  E-2 

-13E-3 

7.4E-4 

2.7E-5 

2.0 

-1.5  E-2 

-8.1E-4 

-13E-5 

1.6E-7 

35 

-12  E-3 

-5.7E-5 

-1.2E-6 

-lOE-8 

5.0 

-5.6  E-4 

-8.4E-6 

-l.OE-7 

-l.lE-9 

2.0 

0.1 

-10  E-1 

3.8E-2 

-7.5E-3 

15E-3 

1.0 

-8.0  E-2 

1.3E-3 

2.7E-4 

-l.OE-5 

2.0 

-11  E-2 

-5.3E-4 

15E-6 

5.4E-7 

35 

-3.9  E-3 

-7.9E-5 

-l.OE-6 

-7.1E-9 

5.0 

-1.1  E-3 

-1.4E-5 

-1.4E-7 

-l.OE-9 

5.0 

0.1 

-1.1  E-1 

1.75E-2 

-3.3E-3 

6.5E-4 

1.0 

-5.1  E-2 

1.0  E-3 

9.8E-5 

-4.0E-6 

2.0 

-1.9  E-2 

-1.3  E-4 

25E-6 

1.7E-7 

35 

-6.2  E-3 

-4.3  E-5 

-1.6E-7 

25E-10 

5.0 

-15  E-3 

-1.5  E-5 

-5.1E-8 

-l.lE-lO 

1.01 

0.1 

-7.3E-1 

^  =  0.707 

4.4E-1 

-14E-1 

1.3E-1 

1.0 

-1.9E-1 

-9.4E-3 

5.9E-3 

4.3E-4 

2.0 

-3.0E-2 

-3.2E-3 

-1.8E-1 

25E-6 

35 

-4.3E-3 

-13E-4 

-9.7E-6 

-3.2E-7 

5.0 

-l.lE-3 

-3.3E-5 

-8.0E-7 

-1.7E-8 
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o  p 

£2(0.  P)^ 

£4(0. 

£4(0,  p)5‘ 

£4(0,  P)4“ 

\  =  0.707  cont. 

2.0  0.1 

-4.0E-1 

15E-1 

-6.0E-2 

2.4E-2 

1.0 

-1.6E-1 

5.1E-3 

2.1E-3 

-1.6E-4 

2.0 

-4.1E-2 

-2.1E-3 

15E-5 

8.7E-6 

35 

-7.9E-3 

-3.2E-4 

-8.2E-6 

-l.lE-7 

5.0 

-2.3E-3 

-5.8E-5 

-l.lE-6 

-1.7E-8 

5.0  0.1 

-2.2E-1 

7.0E-2 

-2.6E-2 

l.OE-2 

1.0 

-l.OE-1 

4.0E-3 

7.8E-4 

-6.4E-5 

2.0 

-3.8E-2 

-5.1E-4 

2.0E-5 

2.7E-6 

35 

-1.2E-2 

-1.7E-4 

-1.2E-6 

3.9E-9 

5.0 

-5.1E-3 

-5.9E5 

-4.1E-7 

-1.8E-9 

If  the  coil  is  to  be  used  as  a  calibration  coil  providing  a  known  uniform  field 
then  an  important  factor  is  the  area  over  whid\  one  requires  field  homogeneity. 
For  example,  if  the  field  homogeneity  is  required  over  a  circular  region  of  1.0  to 
2.0  m  diameter,  then  clearly  a  large  coil  is  required  and  in  order  to  mininrize 
weight  and  costs,  it  is  desirable  to  minimize  a  and  p.  As  an  example,  consider 
a  coil  with  a  radius  of  1.0  m,  coil  thickness  of  1  cm  (a  =  1.01)  and  a  length  of 
20  cm  (P  =  0.1).  (If  a  =  15,  then  the  coil  thickness  is  50  cm  which  leads  to  a 
massive  coil.)  The  field  homogeneity  at  radial  distances  away  from  the  centre 
of  the  coil,  Ji/2)/Ho,  is  as  large  as  7%  at  a  distance  of  r  =  30  cm  and 
deteriorates  to  23%  and  64%  at  50  cm  and  71  cm  away  from  the  coil  centre 
respectively;  see  Table  4.  The  final  coil  dimensions  would  depend  on  the 
required  field  homogeneity  and  acceptable  tolerances. 


Table  4:  Field  homogeneity  for  thin  cross-section  coils  of  internal  radius  1.0  m 


a  =  1.01,  P 

=  0.1  (coil 

thickness,  1  cm;  coil 

length,  20  cm) 

r/a, 

r 

H,  (i  ii/2)/Ho 

H,  (7,  0)/Ho 

0.3 

30  cm 

1.071  (+  7%) 

0.882  (- 12%) 

05 

50  cm 

1.234  (+  23%) 

0.723  (-  28%) 

0.707 

71  cm 

1.637  (+  64%) 

0593  (-  41%) 

a  =  1.01,  P 

=  0.4  (coil 

thickness,  1  cm;  coil 

length,  80  cm) 

r/a^ 

r 

H,  (i  w/2)/Ho) 

H,  (7,  0)/Ho 

0.3 

30  cm 

1.052  (+  5%) 

0.907  (-  9%) 

05 

50  cm 

1.159  (+ 16%) 

0.768  (-  23%) 

0.707 

71  cm 

1.366  (+  37%) 

0.602  (-  40%) 
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2.4  Hart's  Method  for  Distributed  Circular  Currents 


Figure  9  depicts  a  distributed  circular  current  whose  positive  flow  occupies  the 
full  rectangular  area  shown  in  the  shaded  region.  The  coordinates  are 
normalized  to  the  coordinate  P(0,  h)  where  h  is  the  perpendicular  distance  from 
the  2-axis  to  a  point  P.  For  a  distributed  circular  current.  Hart  [2]  obtains  the 
following  expression  for  the  radial  and  axial  magnetic  fields  and  B2 
respectively. 


B 


R 


r  r  E 

R  Z  a 
kJhSg 


RZcos  a(200  A  R  A  Z  A  a) 
[1  +R^+Z^-2RcosaF2 


(23) 


6 


z 


t/k  E  E  E 

R  Z  a 

k}hS2 


(R^  -Rcosa)(200AR  AZ  Aa) 
(1  ♦R2+z2-2RcosaP^ 


(24) 


where  B  is  in  tesla,  distance  is  in  metres,  /  is  the  current  density  in  A  m'^  and 
k  =  10"®  tesla  m  A'^  Sr  and  $2  apply  to  a  uniform  current  distribution  which 
occupies  the  entire  shaded  rectangular  area  of  Figure  9b  from  the  origin 
diagonally  outward  to  point  Q.  The  angular  increment  Aa  was  chosen  to  make 
L  =  Ji/Aa  integral  for  the  semicircle  where  L  is  the  number  of  terms  in  the 
summation  The  summation  is  multiplied  by  2  to  obtain  the  required 
summation  around  the  circle  and  a  is  the  angle  taken  to  the  centre  of  its 
corresponding  increment.  Sr  and  S^  are  numerical  approximations  to  the 
intractable  triple  integrals 

r  r*  (*...drdzda 

Jo  Jo  Jo 


which  have  been  evaluated  digitally  and  presented  in  tabular  form  for 
incremental  values  of  Z  and  R,  see  Hart's  tables  1-A,  B  and  C  [2]. 
Unfortunately,  four-point  interpolation  may  be  necessary  to  deduce  the  Sr  and 
$2  functions  at  the  required  R,  Z  coordinates  from  the  values  given  in  Hart's 
tables. 
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(a)  (b) 


Figure  9:  Distributed  circular  current:  (a)  dimensionless  coordinates, 

(b)  dimensionless  coordinates  normalized  to  h,  i.e.  z  =  Zh,  r  =  Rh,  dz  =  h  dZ  and 
dr  =  h  dR.  The  origin  0  is  taken  at  the  point  at  which  h  meets  the  z-axis. 


Figure  10:  A  distributed  circular  current  of  rectangular  cross  section  in  terms  of 
(a)  dimensional  coordinates,  and  (b)  normalized  coordinates. 
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For  current  distributions  contained  in  the  typical  coil  distribution  shown  in 
Figure  10a,  eqns.  23  and  24  can  be  modified,  (see  Fig.  10b),  resulting  in 


®K,Z  *  ~  ~  *  ^A^R,Z 

where  the  shaded  regions  in  Figure  10b  are  subtracted  from  (contribution  of 
area  DHOF).  However,  the  contribution  from  area  AGOE  has  been  subtracted 
twice,  resulting  in  the  +  term  in  eqn.  25.  Hart's  method  was  used  to 
compute  the  ratio  of  the  axial  field  to  the  central  field  at  various  radial 
distances  from  the  centre  of  a  coil  with  a  =  1.01,  a,  =  100  cm  and  coil  length 
2b  =  40  cm,  100  cm  and  200  cm.  The  results  are  shown  in  Figure  11  and  a 
comparison  of  this  technique  with  Garrett's  series  expansion  approach  is 
summarized  in  Table  5.  The  agreement  is  good  except  for  Hart's  value  at 
h  =  30  cm,  P  =  0.5  which  is  too  high.  The  error  associated  with  this  value 
arises  from  the  two  4-point  interpolations  required  to  extract  the  integral 
from  Hart's  tabulated  values.  Whilst  both  techniques  give  similar  results,  the 
series  expansion  approach  is  much  faster  because  no  numerical  integrations  are 
required  and  the  relatively  simple  E„(a,  P)  error  coefficient  functions  are  readily 
computed.  Accuracy  and  speed  are  also  improved  because  there  is  no  need  to 
resort  to  tables  or  interpolations.  Furthermore,  as  discussed  in  the  following 
section,  the  series  expansion  approach  easily  lends  itself  to  the  computation  of 
fields  from  a  combiiuition  of  coaxial  coils.  Thus  the  series  expansion  approach 
is  far  superior  for  numerical  modelling  studies. 


Table  5:  Comparison  of  Hart's  numerical  integration  procedure  with  Garrett's  series 
expansion  method 


H,  k/D/Hq  for  a  coil  =  100  cm,  Oj  = 
21;  =  40  cm  (P  =  0.2),  100  cm  (P  =  i 

101  cm,  a  =  l.Ol  and  coil  lengths 

05)  and  200  cm  (P  =  1.0) 

P  = 

0.2 

P 

1  =  0.5 

P  = 

1.0 

Garret 

Hart 

Garret 

Hart 

Garret 

Hart 

4  =  0.7 

1.56 

155 

1.28^ 

1.27, 

1.08, 

I.O87* 

fc  =  70  cm 

4  =  0.5 

I.2I5 

I.2I2 

1.13, 

1.13, 

1.04, 

1.04, 

fc  =  50  cm 

4  =  0.3 

I.O67 

1.08j 

1.04, 

1.09, 

I.OI7 

1.02/ 

h  =  30  cm 

*  Obtained  by  linear  interpolation  from  Figure  11 
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H2(h,7r/2)/H 


Figure  11:  Deviation  of  axial  field  HJh,  n/Z)  from  central  field  Hg  computed  from 
Hart's  tables  of  functions  for  coil  length  40  cm,  100  cm  and  200  cm. 
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3.  Fields  for  Two  Coaxially  Displaced  Solenoids 

3.1  Introduction 


Using  the  superposition  principle  and  symmetry  arguments,  the  axial  field 
along  the  coil  axis,  H,(z,  0),  outside  or  within  a  coil  or  between  two  or  more 
coils  separated  by  a  gap  can  always  be  expressed  as  a  sum  of  central  fields. 

For  example,  the  end  field  0)  of  any  coil  is  half  the  central  field  of  a 

coil  twice  the  length.  The  field  Hjly,  0)  for  y  >  p,  i.e.  external  to  the  coil  can  be 
computed  from  the  central  fields  of  two  sub-coils  into  which  the  original  coil  is 
partitioned,  see  Figure  12.  Normalized  distances  are  used  for  the  sub-coils. 
Sub<oil  A  has  geometry  (a,  P  +  y)  and  sub-coil  B  has  geometry  (a,  y  -  p).  Thus 


H,(y,0) 


(26) 


where  from  eqn.  11,  the  central  fields  of  sub-coils  A  and  B  are  given  as 

Ho"*  -  ^Xa,F(a,p  +y) 


and  thus 


where  Hq  is  the  central  field  Hj(0,0)  of  the  original  coil  with  geometry  (a,  P),  i.e. 

P).  If  the  axial  field  is  required  at  a  point  along  the  axis  which  lies 
within  the  coil  boundary,  then  eqn.  10  can  be  used  to  modify  the  term  F(a,  y  - 
p)  appearing  in  eqn.  28,  see  Figure  13.  Thus 

mH  f(tt»P +y)  ■•■F(o,p -y)  129) 

lH,(y,0)J^,p  - 

From  eqn.  8,  eqns.  28  and  29  reduce  to  eqn.  19. 
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Figure  12:  Partitioning  of  main  coil  into  sub  coQ  A  and  sub  coil  B  allows  the  external 
field  HJz,  0)  to  be  evaluated  from  the  central  field  of  sub  coils  A  and  B.  y  =  z/uj. 


3.2  Axial  Fields  from  Multiple  Coils 


The  previous  analysis  by  s)aninetiy  and  superposition  arguments  can  also  be 
applied  to  more  than  one  coil,  for  example,  a  pair  of  coils  [1].  The 
homogeneity  between  a  pair  of  coils  can  then  be  analysed  by  the  Legendre 
series  expansion  techixique. 
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Hz(z/a^,0)=  Ho/2  ♦Ho/2 


Figure  13:  Partitioning  of  main  coil  into  sub  coU  A  and  sub  coil  B  allows  the  internal 
field  Hj[z,  0)  to  be  evaluated  from  the  central  field  of  sub  coils  A  and  B.  y  =  z/aj. 


For  the  pair  of  coils  described  in  Figure  14,  the  axial  field  along  the  axis  is 


H,(y,0) 


*  [F(a,2p+P,+Y)-F(a,p^+T) 

*F(a,2p+p^-Y)-F(a,P^-Y)] 


(30) 


where  P^  =  g/2ay  F(a,  C)  is  defined  by  eqn.  8  and  H^{-  y,  0)  =  HJiy,  0).  The 
field  at  the  centre  of  the  gap  can  then  be  expressed  as 
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H,(0,0)  -  j^Xa^  [F(a2p+pp-f  (a,pp] 

+  (a,p,pp 


(31) 


Equation  31  is  equivalent  to  a  coil  of  length  4p  +  2p^  from  which  a  coil  of  length 
2P2  has  been  subtracted. 


Q/a-j-f  z/ai 


2^-t'  G/a-j  ♦  2/a.| 


-  Fla,  g/2ai  +  z/ai) 


+  F (a,  2^+  9/23-1+2/31) 


G  =  g/2 
/5g  =  g/2ai 


♦  F|a2/3  +  g/2ai-z/3i) 


-Flag  /2ai  -  z  /ai) 


Figure  14:  Partitioning  of  a  pair  of  coaxial  coils  into  four  sub  coils,  y  =  z/aj. 
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The  deviation  of  the  axial  field  from  the  central  field,  HJiy,  0)/HJi0J0),  of  a 
split  coil  pair  as  one  moves  away  from  the  centre  of  the  gap  can  be  computed 
from  eqn.  30.  The  results  for  a  fat  coil  (a  =  3,  P  =  1.0)  and  a  thin  filamentaiy 
loop  coil  (a  =  1.01,  P  =  0.1)  are  shown  in  Figures  15, 16  and  17,  where  the  gap 
width  parameter  is  varied  from  zero  to  2.0  (gap  width  fotu*  times  the  internal 
radius).  Rgures  15  and  16  show  two  different  coils  with  the  same  settings. 
Bgure  17  focuses  on  the  homogeneous  region  of  Figure  16.  In  Figure  15,  there 
are  certain  gap  settings  (p  -  0.25)  where  the  field  homogeneity  along  the  axis  is 
within  a  few  percent  of  Hg  out  to  z/a,  =  1.0.  A  similar  result  applies  to  the 
"filamentary"  loop  coil,  except  that  the  region  of  field  homogeneity  holds  for  a 
different  value  and  does  not  extend  out  to  z/a^  =  1.0  as  for  the  fat  coil.  For 
example,  for  the  thin  foil.  Figure  17  reveals  that  the  optimum  field  homogeneity 
is  obtained  at  about  =  0.45  up  to  z/fl|  =  0.4  (HJiy,  0)/Hq  =  0.9982).  If  a  = 

1.01,  P  =  0.1  and  =  100  cm,  then  the  optimum  gap  setting  is  approximately  90 
cm  and  z  =  40  cm  for  field  homogeneity  to  within  0.2%  of  the  central  field. 
Figures  15  and  16  also  show  the  axial  field  deviation  for  zero  gap  width.  Here 
the  field  drops  off  monotonically  away  from  the  coil  centre. 

The  axial  field  within  a  region  of  space  around  the  centre  of  a  single  coil  is 
less  imiform  along  the  z  axis  (6  =  0)  than  radially  outwards  (0  =  n/2),  see 
section  23  and  Table  4.  If  the  same  conclusion  applied  to  a  pair  of  coils,  then 
eqn.  30  and  Figures  15  to  17  can  be  used  to  quantitatively  set  the  upper  limits 
of  field  homogeneity  within  a  spherical  region  about  the  centre  of  a  split  coil 
pair.  However  this  is  not  always  the  case  and  depends  on  the  gap  separation 
and  size  of  coils.  Choosing  the  appropriate  P^  for  good  axial  field  uniformity 
along  the  z  axis  (P^  =  0.45,  Fig.  17),  eqns.  30  and  31  are  evaluated  and  the 
uniformity  H,(z/flj,  0)/Hq  for  certain  z/a^  is  shown  in  Figure  18.  (Here  the 
values  in^cate  the  percentage  deviation  from  Hg  at  r/a^  and  z/Oi  =  0.3,  0.4,  0.5, 
0.6,  0.8.  The  deviations  in  brackets  are  those  obtained  from  the  Legendre 
expansion  technique,  see  section  3.3).  Note  that  the  individual  coils  (a  =  1.01, 

P  =  0.1)  are  the  same  as  the  example  given  in  Table  4  and  by  comparison,  the 
split  coil  pair  offers  better  field  imiformity  by  more  than  an  order  of  magnitude 
up  to  z/ui  =  0.5.  For  example  at  z/a^  =  0.5,  the  deviation  from  the  central  field 
is  28%  for  a  single  coil  compared  to  2%  for  the  split  coil  pair  and  from  Figure  7 
at  a  =  1.01,  a  P  value  of  1.4  is  required  to  achieve  2%  single  coil  uniformity. 

This  implies  a  coil  length  of  2.8  m  as  opposed  to  0.2  m.  There  are  several  other 
interesting  features  revealed  in  Figiue  18.  The  Legendre  expansion  technique 
which  includes  terms  up  to  eighth  order  yields  results  in  excellent  agreement 
with  the  superposition  expressions  (eqns.  30  and  31).  The  discrepancy  at  larger 
z/Hi  arises  from  the  contributions  of  higher  order  terms.  For  example,  the 
eighth-order  term  accormts  for  10%  of  the  computed  HJiy  =  0.8,  0)/Hg  value 
(0.86480).  Another  feature  is  that  the  field  uniformity  in  the  6  =  n/2  plane  is 
also  improved  by  using  a  pair  of  coils  rather  than  a  single  coil  (see  Table  4). 

The  field  uniformity  for  the  coil  pair  is  better  in  the  z  direction  than  the  0  =  n/2 
plane.  Similarly,  for  a  =  3.0,  P  =  1.0  (Fig.  15),  the  same  trend  obtains  at  P^  = 

0.2  but  reverses  at  p^  =  1.0. 
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Figure  15:  Deviation  of  axial  field  H/zlaj,  0)  from  field  at  gap  centre  for  a  coil  pair 
(a  =  3.0,  p  =  1.0). 
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H2(7  /ai,o)/H 


Figure  16:  Deviation  of  axial  field  H/zlaj,  0)  from  field  at  gap  centre  for  a  coil  pair 
(a  =  1.01,  p  =  0.1). 
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z/ai 


Figure  17;  Deviation  of  axial  field  H/z/oj,  0)  from  field  at  gap  centre  for  a  coil  pair 
(a  =  1.01,  p  =  0.1). 
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a  =  101 
fi  =  0-1 

/3g=  0-^5 


Figure  18:  Regions  offidd  un^ormity  Hjlp/ap  Q  =  0,n/2  fm-  a  split  coil 
a  =  1.01,  p  =  0.1,  p,  =  0.45. 

3.3  Fields  Within  a  Central  Region  for  a  Pair  of  Coils 

The  Legendre  expansion  technique  is  combined  with  the  superposition  principle 
for  multiple  coils  in  order  to  obtain  inhomogeneities  in  the  central  region  of  the 
gap  formed  by  a  pair  of  coaxially  displaced  coils.  Thus  the  expressions  for  the 
axial  fields  HJly,  0)  and  H,(^,  n/2),  eqns.  20  and  21,  will  contain  for  each  term  a 
corresponding  term  representing  the  compensating  coil  contribution  which  for 
the  pair  of  coils  is  the  imaginary  gap  or  void  coil.  There  is  a  degree  of 
freedom  in  choosing  the  location,  inner  and  outer  radius  and  relative  current 
density  of  the  compensation  coil.  However,  the  number  of  terms  in  the  series 
expat\sion  of  the  field  generated  by  the  main  coU  that  can  be  cancelled  by  the 
terms  arising  from  the  compensating  coil  is  determined  by  the  number  of 
independent  variables  of  the  comfjensating  coil.  Thus  for  the  separated  coil 
pair  in  Figure  14,  the  superimposed  "negative"  coil  produces  the  gap  (upper 
and  lower  two  sub-coils  in  Fig.  14)  and  since  the  main  coil  and  compensating 
coil  have  the  same  o,  the  only  free  variable  is  p^,  the  gap  separation. 
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Consequently,  only  the  £2(0,  P)  term  in  the  series  expansion  may  be  cancelled 
leading  to  a  fourtti-order  compensated  coU.  On  the  other  hand  if  the  main  coil 
is  configured  to  have  a  notch  either  on  the  inside  or  outside  of  the  coil  winding, 
positioned  symmetrically  about  the  mid-plane,  then  both  (s  the  depth  of  die 
notch)  and  (a  the  length  of  the  notch)  may  be  varied.  Consequendy,  both 
the  £2(01,  P)  and  £4(0,  P)  terms  may  be  cancelled  by  choosing  the  appropriate 
(a^  Pp  in  order  to  achieve  even  higher  field  homogeneities.  This  configuration 
leads  to  a  sixth-order  compensated  coil  [6]. 

Let  a,  P  represent  the  geometry  parameters  for  the  main  coil  including  the 
gap  as  shown  in  Appendix  A,  (unlike  Hg.  14)  and  let  P^  refer  to  the 
compensating  coil  (gap).  Then,  from  eqns.  20  and  21, 


HPy,0)  -  (f  (c^p)  -fpa^pp 

+  [Fo£2(a,p)-f^£2(a^ppiT^ 

+  [fo£4(a,p)-f^£4(a^pp]Y* 

*[Fo£pa,p)-f^£g(a^,pp]V» 

+  ...), 


Hptn/2)  =  jM 

-.l[Fo£2(a,p)-F^£2(a^Pp]^2 

♦  |(Fo£4((x,p)-F^£,(a^pp]5< 

-^[Fo£4(a,p)-F^£,(aj,ppie 

>...) 


where  Fj£„((Xi,  pj)  is  the  product  of  the  field  factor  (eqn.  8)  and  eqn.  18  (see 
Table  2)  both  evaluated  at  Oj,  P,.  If  Fo£2(a,  P)  =  F^2^a^  Pp,  the  main  coil  will 
be  compensated  to  fourth  order.  This  defines  a  Helmholtz  coil  pair. 

A  graphical  and  hence  approximate  solution  can  be  obtained  by  plotting 
F£2(a,  P)  as  a  function  of  p  for  various  values  of  a.  The  results  are  shown  in 
Figures  19  to  22  which  can  be  readily  used  for  determining  the  correct  gap 
widths.  The  expanded  left  hand  portion  of  Figures  19  and  21  appear  in 
Figures  20  and  22  respectively  and  are  required  for  determining  the  correct  P^ 
value  for  fourth  order  compensation.  For  example,  for  a  =  3,  P  =  2, 

FE2(3,  2)  =  -  0.153  ■»  ££2(3, 0.193);  thus  a  gap  setting  P^  =  0.193  will  yield  the 
maximum  central  homogeneity.  It  is  evident  from  Figures  19  to  22  that  for  a 
given  a,  the  further  away  that  P  is  from  the  minimum  of  the  ££2(0,  P)  curve, 
then  the  smaller  the  required  P^  value  for  fourth-order  compensation.  This 
leads  to  fat  coils  with  a  narrow  spacing  (P^  «  P)  coils.  On  the  other  hand,  if  P 


37 


is  chosen  close  to  the  minimum  then  is  a  large  proportion  of  P  and  this 
configuration  gives  rise  to  thin  pancake-type  coOs  which  are  lighter  and  more 
economical  than  the  fat  «  p)  coils.  As  an  example,  let  a  =  1.01  and 
a-i  =  100  cm.  When  P  =  1.0,  fE2{a,  P)  =  -  6.63  x  10^  and  from  Figure  22,  p^  is 
020,  ^£2(0,  pp  =  -  6.65  X  10^;  this  leads  to  the  coil  configuration  shown  in 
Figure  23a.  If,  however,  p  =  0.65,  ££2(0,  Pp  =  -  1.00  x  10'^  then  p^  =  0.39  and 
one  obtains  the  configuration  described  in  Figure  23b. 


See  Fig  21 


Figure  19:  FEjia,  as  a  function  of  ^  for  a  ^  1.01, 1.05, 1.1, 1.2, 13,  1.4, 1.5, 1.7, 
2.0,  23, 3.0,  4.0  and  5.0. 
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Figure  20:  Expanded  portion  of  Figure  19. 
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Figure  21:  FE^^a,  as  a  function  of  ^  for  a  =  1.01,  1.02, 1.03,  1.04  and  1.05. 
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Figure  23:  Fourth  order  compensated  coib  (scale  1:32). 


Another  interesting  feature  displayed  by  Figures  19  and  21  is  that  for  small 
values  of  a  12)  the  minimum  F£2(a,  P)  value  occurs  at  ^  0.5.  Correct 
compensation  at  the  ££2(0,  P)  minimum  implies  «  p.  llus  leads  to  the 
approximate  result  that  for  Helmholtz  coils,  the  gap  width  is  equal  to  the 
internal  radius  of  short  coils.  As  a  increases,  the  ££2(0,  p)  minimum  occurs  at 
larger  P  values.  At  a  =  5.0,  the  ££2(0,  P)  minimum  obtains  at  P  =  0.8  (P^  »  P) 
and  the  required  gap  width  for  Helmholtz  coils  is  1.6a,.  Figures  19  and  21  are 
also  useful  for  determining  the  minimiun  values  of  P  for  a  given  a  to  allow  for 
fourth-order  compensation.  If  P  S  0.8  for  a  coil  having  a  =  5.0,  then  fourth- 
order  compensation  is  not  possible. 
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3.4  Numerical  Computation  of  Compensated  Fourth 
Order  Coils 

In  the  previous  section,  computer  geirerated  plots  were  used  to  provide 
graphical  solutions  for  choosing  tlve  coil  parameters  P  and  for  a  given  a. 
However,  it  is  also  useful  to  automatically  compute  the  correct  coil  parameters 
for  obtaining  compensated  coils.  This  allows  a  more  comprehensive  study  to 
be  undertaken  and  is  also  useful  for  tabulating  error  limit  tables  (see  section 
3.5). 


3.4.1  The  Newton-Raphson  Method  for  Determining 

To  compute  HJiy,  0)  and  HJi%,  n/2)  (eqns.  32,  33)  to  fourth-order,  it  is  necessary 
to  solve  (see  Table  2) 


1 

2F 


/ 

\  1 

f 

\ 

1 

\/2 

1 

1 

2 

i/2 

(1  +32)3/2 

0^  +  3^. 

j 

’W, 

(1  +3^^^^ 

^  * 

“3  ^ 

> 

(34) 


for  3^  in  terms  of  a  and  3  (a  =  for  Helmholtz  coils).  Equation  34  cannot  be 
solv^  analjftically  and  so  must  be  solved  numerically  using,  for  example,  the 
Newton  Raphson  technique  for  the  solution  of  nonlinear  equations  [7].  This 
technique  is  well-suited  to  this  case  because  the  derivative  of  the  function  can 
be  expressed  analytically.  However,  because  of  the  geometry  of  the  method, 
the  Newton-Raphson  method  can  be  erratic  in  regions  where  the  function  has  a 
small  slope. 

The  Newton-Raphson  method  attempts  the  solution  of  the  equation  f  (x)  =  0 
where  if  x,  is  an  approximation  to  x,  then  a  better  approximation  is 


or 


(35) 


AX;F'(Xi) +F(x,.)  =0 


(36) 


where  F'(x,)  is  the  first  derivative  of  F(x,),  in  this  case,  3)  -  Pg)- 

If  AXj  (=  x,+i  -  Xj)  is  less  than  a  given  tolerance,  i.e.  Ax,  approaches  zero,  then  X; 
is  an  approximate  solution. 
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3.5  Error  Limits  for  Helmholtz  Coils 


As  pointed  out  previously,  the  terms  in  eqns.  16  and  17  have  their  largest  value 
whm  6  =  0,  P„(l)  =  1,  ai^  the  maximum  value  of  the  terms  for  the  axial  field  is 
given  by  eqn.  20  or  32.  When  the  second  term  in  eqn.  32  is  zero,  there  lemams 
only  the  fourth  order  term  and  the  deviation  from  the  central  field  H^/Hg  is 


H,(y,0)  -  1  + 


FoE^(a,p)-f  E,(a,p 


Fg(a,p)-fJa,pj 


l+A-/ 


+  ... 


(37) 


where 


(38) 


and  if  A  is  the  percentage  change  in  HJiy,  O)  from  Hg  then 


z 


(39) 


determines  how  far  from  the  origin  one  can  move  before  the  field  deviates  by 
more  than  a  given  percentage  A.  This  assumes  that  sixth-order  and  higher 
contributions  are  negligible.  The  evaluation  of  z/a^  in  eqn.  39  is  readily 
determined  using  the  Newton-Raphson  method  to  automatically  compute  the 
required  P^  for  cancellation  of  the  second  order  term.  Some  sample 
computations  for  "short"  and  "long"  coU  pairs  are  presented  in  Table  6.  Thus, 
depending  upon  the  acceptable  errors,  a  significant  portion  of  the  coil  may  be 
used  to  provide  an  approximately  uniform  magnetic  field  for  calibration 
purposes. 
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Table  6:  Axial  field  error  limits  for  fourth-order  coils 


2/«X- 


a 

P 

P* 

r/cm 

0.1% 

1.0% 

5.0% 

1.01 

1.0 

0.197 

200 

39.4 

0201 

0.357 

0334 

0.6 

0.415 

120 

83 

0.174 

0309 

0.463 

1.03 

1.0 

0.202 

200 

40.4 

0.202 

0359 

0337 

03 

0.424 

120 

84.8 

0.173 

0312 

0.467 

1.05 

1.0 

0.207 

200 

41.4 

0.203 

0361 

0340 

0.6 

0.433 

120 

86.6 

0.173 

0315 

0.471 

1.1 

1.0 

0.221 

200 

44.2 

0.206 

0366 

0347 

0.6 

0.454 

120 

90.8 

0.181 

0322 

0.481 

13 

1.2 

0.236 

240 

472 

0241 

0.428 

0.640 

0.8 

0.433 

160 

86.6 

0212 

0377 

0364 

2.0 

1.3 

0279 

260 

55.8 

0265 

0.472 

0.706 

0.8 

0343 

160 

108.6 

0234 

0.416 

0.623 

4.0 

13 

0.363 

300 

72j6 

0.321 

0371 

0.853 

0.9 

0.666 

180 

133.2 

0288 

0313 

0.767 

*  I  and  tg  refer  to  overall  coil  length  including  gap,  and  gap  width 
respectively  for  an  inner  radius  of  100  cn\. 

**  z/uj  determines  how  far  from  the  origin  one  can  move  before  the  field 
deviates  by  more  than  the  given  percentage. 


3.6  Axial  Fields  for  Uncompensated  Coil  Pairs 


As  shown  in  Figures  15  to  17,  it  is  possible  to  forgo  maximum  field 
homogeneity  for  the  sake  of  extending  z/a-i  such  that  the  deviation  from  the 
central  field  is  within  certain  error  limits.  For  example,  the  =  0.45  curve  in 
Figure  17  extends  further  out  than  the  P^  =  0.4  curve,  etc.  It  is  thus  desirable 
to  evaluate  0)  and  Hj(^,  n/2)  for  a  coil  pair  when  the  second-order  term 
has  not  been  cancelled  out,  i.e.  evaluate  eqns.  32  and  33.  A  Fortran  source 
code  listing  for  the  evaluation  of  these  equations  is  described  briefly  in 
Appendix  A;  the  programme  also  generates  two  plot  files  for  graphing 
HJiy,  0)/Hq  and  H/.^,  n/0)/Hg  and  a  scratch  file  listing  the  second-,  fourth-, 
sixth-  and  eighth-order  contributions  so  that  the  convergence  of  eqns.  32  and  33 
can  be  monitored.  Thus  one  can  readily  examine  the  deviation  of  the  axial 
field  from  the  central  field,  both  radially  and  along  the  z-axis,  for  a  range  of  P^ 
values,  including  the  one  required  for  fourth-order  compensation. 
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3.6.1  Experimental  Design  Example 


As  an  example,  consider  a  split  coil  with  two  adjustable  lengths  where  a  =  1.01 
and  P  =  030  and  0.55.  The  confuted  axial  field  deviations  HJizfay  0)/Hq  and 

H, (r/fli,  k/D/Hq  for  various  gap  settings  are  shown  in  Figures  24  and  25.  As 
is  evident  from  these  graphs,  as  approaches  05,  the  correct  fourth-order 
compensation  coU  setting,  better  field  homogeneity  is  obtained.  With  regards 
to  HJr/a^,  kIIUHq,  it  is  also  dear  from  these  figures  that  other  values  of  P^  still 
give  a  very  small  deviation  (albeit  a  larger  one)  from  the  central  field  at  far 
greater  displacements  from  the  gap  centre.  There  is  also  very  little  variation 
between  the  results  shown  in  Figures  24  and  25  and  those  obtained  from  a  = 

I. 0(^,  1.02  and  1.03  which  suggests  that  it  is  advantageous  to  use  coils  with 
small  a  to  allow  for  easier  winding  and  less  overall  weight.  The  use  of  low  a 
values  becomes  important  when  field  homogeneity  is  required  over  a  large 
area. 

In  Figure  24,  the  axial  field  Hjirla^,  n/2)  at  P^  =  0.35  does  not  deviate  by  more 
than  1%  from  Hq  for  displacements  of  up  to  60%  of  the  internal  radius  of  the 
coils.  The  axial  field  0)  is  also  uniform  to  within  1%  of  Hq  for 

displacements  of  up  to  0.15ai  each  way  along  the  z-axis.  Thus  for  a  standard 
internal  radius  of  100  cm  the  field  is  expected  to  be  uniform  to  within  1%  of  Hq 
for  a  region  of  space  bounded  by  a  disc  (sensor  coil)  having  a  radius  of  60  cm 
and  a  tMckness  of  30  cm.  This  was  verified  experimentally.  If  better  field 
uniformity  is  required  over  a  smaller  region  of  space,  then  one  can  increase  the 
size  of  the  gap  using  the  original  peur  of  coils  (see  Fig.  25)  where  P  has  now 
been  increased  from  05  to  0.55  and  P^  =  0.40.  Now,  HJir/uy  n/2)  [Hj(z/aj,  0)] 
is  uniform  to  within  0.2%  of  Hg  for  displacements  up  to  r/uj  =  0.4  [z/aj  =  0.1]. 
Tables  7  and  8  contain  the  information  presented  in  Figures  24  and  25  in  tabular 
form  to  allow  for  better  quantitative  estimates  of  field  uniformity  as  these  coil 
coirfigurations  will  be  used  experimentally. 

The  accuracy  of  the  uniformity  can  be  checked  from  the  convergence  of  the 
series  expansion  given  by  eqns.  32  and  33.  The  even  error  coefficients 
FoE„(a,  P)  -  FgE„(a,  Pp  up  to  eighth  order  for  the  coaxially  displaced  coil  pairs 
given  in  Figures  24  and  25  are  presented  in  Table  9  for  p^  =  0.35,  P  =  05  and  P^ 
=  0.4,  P  =  055.  For  a  =  1.01,  P  =  055,  P^  =  0.4,  the  results  in  Table  9  show  that 
up  to  r/uj  =  0.4  the  series  has  converged  rapidly  and  that  it  is  not  necessary  to 
go  beyond  the  fourth-order  term.  On  the  other  hand,  for  a  =  1.01,  P  =  05,  P^  = 
035  it  is  necessary  to  go  beyond  the  fourth-order  term  to  obtain  conveigence 
for  r/«i  S  0.65,  z/a,  S  0.65.  From  Table  7,  H,(0.65,  n/2)  =  0.97  Hg  and  from 
eqn.  33  and  Table  9,  the  eighth-order  contribution  is  0.98%  of  the  total  Hj(0.65, 
Jt/2/Ho.  Since  higher  order  terms  contribute  less,  the  uniformity  of  HJr/a-i, 
n/2)/Ho  to  better  than  3%,  for  r/aj  up  to  0.65,  is  not  affected  by  truncating  the 
series  at  the  eighth-order  term. 
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(2/ai.o)/Ho  Hz(r/ai.7r/2)/H 


Figure  24:  Deviation  of  axial  fields  HJirlUi,  nfl)  and  HJlzIaj,  0)  from  the  central  field 
Hq  for  a  coil  pair:  a  =  1.01,  P  =  0.50  and  P^  =  0.48,  0.45, 0.40,  0.35,  03  and  0.22. 
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Table  7:  Axial  fidd  deviations  from  the  central  field  at  various  gap  settings  for 
a  =  1.01,  p  =  050 


p,=0.48 

II 

o 

pj  =  0.40 

P,  =  035 

P,  =  030 

II 

o 

r/oj 

l.OOOlE+00 

l.OOOlE+00 

1.0OO3E-hOO 

1.000^-h00 

1.0005E-fOO 

1.0007E+00 

0.05 

1.0002E-h00 

1.0005E+00 

1.0010E-»00 

1.0015E+00 

1.0021E+00 

1.0029E-fO0 

0.10 

1.0(X)3E+00 

1.00l0E-h00 

1.0022E-f00 

1.O034B-HOO 

1.0046E+00 

1.0065E+00 

0.15 

1.0003E+00 

1.0015E+00 

1.0036E+00 

1.0058E+00 

1.0080E-h00 

1.0115E+00 

030 

9.9979E-01 

1.001 7E+00 

1.0050E+00 

1.0085E+00 

1.0121E+00 

1.0177E+00 

035 

9.9852E-01 

1.0013E+00 

1.0062E+00 

1.0113E+00 

1.0166E+00 

1.0251E+00 

030 

9.9599E-01 

9.9985E-01 

1.0066E-e00 

1.0138E-i^ 

1.0212E+00 

1.0333E+00 

035 

9.9158E-01 

9.9668E-01 

1.0057E+00 

1.0154E-hOO 

1.0255E+00 

1.0420E-fOO 

0.40 

9.8444E-01 

9.9097E-01 

1.0026E-f00 

1.0152E+00 

1.0286E+00 

1.0508E+00 

0.45 

9.7352E-01 

9.8164E-01 

9.9631E-01 

1.0123E+00 

1.0295E+00 

1.0587E+00 

050 

95745E-01 

9.6726E-01 

9.8524E-01 

1.0052E+00 

1.0270E+00 

1.0647E+00 

055 

9.3446E-01 

9.4602E-01 

9.6753E-01 

9.9192E-01 

1.0191E+00 

1.0674E+00 

0.60 

9.0237E-01 

9.1558E-01 

9.4069E-01 

9.6991E-01 

1.0033E+00 

1.0644E+00 

0.65 

85837E-01 

8.7297E-01 

9.0147E-01 

9.3577E-01 

9.7624E-01 

1.0530E+00 

0.70 

7.9901E-01 

8.1443E-01 

8.4573E-01 

8.8511E-01 

9.3345E-01 

1.0292E+00 

0.75 

7.1997E-01 

7.3526E-01 

7.6816E-01 

8.1221E-01 

8.6912E-01 

9.8770E-01 

0.80 

6.1591E-01 

6.2958E-01 

6.6207E-01 

7.0979E-01 

7.7574E-01 

93180E-01 

0.85 

4.8031E-01 

4.9010E-01 

5.1903E-01 

5.6857E-01 

6.4365E-01 

832S7E-01 

0.90 

3.0518E-01 

3.0786E-01 

3.2855E-01 

3.7688E-01 

4.6063E-01 

6.7863E-01 

0.95 

8.0824E-02 

7.1856E-02 

7.7606E-02 

1.2014E-01 

2.1128E-01 

4.7554E-01 

1.00 

H,(z/a„  0)/Ho 

P,=0.48 

p,  =  0.45 

Pj  =  0.40 

II 

o 

II 

o 

II 

o 

2/«l 

9.9987E-01 

9.9973E-01 

9.9947E-01 

9.9921E-01 

9.9895E-01 

9.9853E-01 

0.05 

9.9941E-01 

9.9883E-01 

9.9782E-01 

9.%78E-01 

9.9574E-01 

9.9410E-01 

0.10 

9.9838E-01 

9.9708E-01 

9.9484E-01 

9.9255E-01 

9.9024E-01 

9.8664E-01 

0.15 

9.9641E-01 

9.9413E-01 

9.9023E-01 

9.8624E-01 

9.8226E-01 

9.7606E-01 

030 

9.9303E-01 

9.8954E-01 

9.8359E-01 

9.7756E-01 

9.7154E-01 

9.6226E-01 

035 

9.8772E-01 

9.8282E-01 

9.7453E-01 

9.6616E-01 

95787E-01 

9.4516E-01 

030 

9.7993E-01 

9.7349E-01 

9.6263E-01 

95176E-01 

9.4105E-01 

9.2475E-01 

035 

9.6920E-01 

9.6112E-01 

9.4757E-01 

9.3410E-01 

93093E-01 

9.0105E-01 

0.40 

95516E-01 

9.4540E-01 

93912E-01 

9.1305E-01 

8.9746E-01 

8.7413E-01 

0.45 

93760E-01 

9.2615E-01 

9.0716E-01 

8.8854E-01 

8.7060E-01 

8.4406E-01 

050 

9.1653E-01 

9.0339E-01 

8.8169E-01 

8.6052E-01 

8.4028E-01 

8.1075E-01 

055 

8.9218E-01 

8.7730E-01 

85274E-01 

8.2889E-01 

8.0625E-01 

7.7386E-01 

0.60 

8.6502E-01 

8.4819E-01 

83033E-01 

7.9329E-01 

7.6785E-01 

7.3245E-01 

0.65 

8.3573E-01 

8.1642E-01 

7.8423E-01 

75289E-01 

7.2367E-01 

6.8465E-01 

0.70 

8.0511E-01 

7.8228E-01 

7.4372E-01 

7.0595E-01 

6.7110E-01 

63713E-01 

0.75 

7.7398E-01 

7.4576E-01 

6.9727E-01 

6.4937E-01 

6.0574E-01 

55446E-01 

0.80 

7.4298E-01 

7.0626E-01 

6.4201E-01 

5.7805E-01 

53057E-01 

45820E-01 

0.85 

7.1235E-01 

6.6220E-01 

5.7313E-01 

4.8399E-01 

4.0497E-01 

3.2587E-01 

0.90 

6.8153E-01 

6.1054E-01 

4.8312E-01 

35532E-01 

2.4350E-01 

1.3%2E-01 

0.95 

6.4881E-01 

5.4614E-01 

3.6078E-01 

1.7499E-01 

1.4403E-02 

-1.2537E-01 

1.00 
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Table  8:  Axial  field  deviations  from  the  central  field  at  various  gap  settings  for 
a  =  1.01,  p  =  055 


p,  =050 

p,  =  0.45 

p,  =  0.40 

P,  =  055 

P,  =030 

P,  =  032 

r/«i 

9.9989E-01 

1.0000E-fO0 

1.0001E-h00 

1.0003E+00 

1.0OO4E-hOO 

1.0006E-i^ 

0.05 

9.9954E-01 

l.OOOOE+00 

1.0005E+00 

l.OOllE+00 

1.0O16E-HOO 

1.0025E+00 

0.10 

9.9885E-01 

9.9995E-01 

l.OOllE+00 

1.0023E+00 

1.0036E+00 

1.0055E+00 

0.15 

9.9762E-01 

9.9961E-01 

1.001 7E+00 

1.0039E+00 

1.0061E+00 

1.0097E+00 

030 

9.9562E-01 

9.9875E-01 

1.0021E+00 

1.0056E+00 

1.0091E+00 

1.0149E+00 

035 

9.9247E-01 

9.9702E-01 

1.0019E+00 

1.0070E+00 

1.0123E+00 

1.0208E+00 

030 

9.8770E-01 

9.93%E-01 

1.0007E+00 

1.0078E+00 

1.0153E+00 

1.0274E+00 

035 

9.8069E-01 

9.8892E-01 

9.9786E-01 

1.0074E+00 

1.0175E+00 

1.0341E+00 

0.40 

9.7065E-01 

9.8111E-01 

9.9260E-01 

1.0050E+00 

1.0182E+00 

1.0403E+00 

0.45 

93659E-01 

9.6948E-01 

9.8383E-01 

9.9956E-01 

1.0165E+00 

1.0454E+00 

030 

9.3725E-01 

95270E-01 

9.7016E-01 

9.8964E-01 

l.OllOE+00 

1.0481E+00 

035 

9.1112E-01 

92909E-01 

9.4981E-01 

9.7341E-01 

9.9983E-01 

1.0470E+00 

0.60 

8.7630E-01 

8.9652E-01 

9.2046E-01 

9.4846E-01 

9.8064E-01 

1.0399E+00 

0.65 

8.3048E-01 

85236E-01 

8.7917E-01 

9.1165E-01 

95022E-01 

1.0240E+00 

0.70 

7.7088E-01 

7.9330E-01 

8.2221E-01 

85893E-01 

9.0443E-01 

9.9554E-01 

0.75 

6.9412E-01 

7.1530E-01 

7.4490E-01 

7.8515E-01 

8.3788E-01 

9.4950E-01 

0.80 

5.9617E-01 

6.1337E-01 

6.4137E-01 

6.8376E-01 

7.4368E-01 

8.7942E-01 

0.85 

4.7221E-01 

4.8142E-01 

5.0432E-01 

5.4649E-01 

6.1304E-01 

7.7698E-01 

0.90 

3.1652E-01 

3.1206E-01 

3.2477E-01 

3.6305E-01 

4.3490E-01 

6.3151E-01 

0.95 

1.2235E-01 

9.6354E-02 

9.1663E-02 

1.2064E-01 

1.9540E-01 

4.2958E-01 

1.00 

H,(2/a„  0)/Ho 


II 

o 

P,=0.45 

II 

o 

o 

II 

p,  =030 

P,  =  032 

z/fli 

1.0002E+00 

9.99%E-01 

9.9970E-01 

9.9944E-0* 

9.991 7E-01 

9.9875E-01 

0.05 

1.0007E+00 

9.9974E-01 

9.9873E-01 

9.9768E-01 

9.9663E-01 

9.9496E-01 

0.10 

1.0013E+00 

9.9912E-01 

9.9686E-01 

9.9455E-01 

9.9221E-01 

9.8855E-01 

0.15 

1.0015E+00 

9.9770E-01 

9.9377E-01 

9.8974E-01 

9.8569E-01 

9.7937E-01 

0.20 

1.0008E+00 

9.9503E-01 

9.8901E-01 

9.8289E-01 

9.7678E-01 

9.6728E-01 

0.25 

9.9873E-01 

9.9054E-01 

9.8213E-01 

9.7362E-01 

9.6517E-01 

93213E-01 

0.30 

9.9453E-01 

9.8370E-01 

9.7625E-01 

9.6156E-01 

93060E-01 

9.3383E-01 

0.35 

9.8764E-01 

9.7399E-01 

9.6016E-01 

9.4637E-01 

9.3284E-01 

9.1232E-01 

0.40 

9.7760E-01 

9.6102E-01 

9.4435E-01 

93784E-01 

9.1177E-01 

8.8760E-01 

0.45 

9.6409E-01 

9.4456E-01 

9.2505E-01 

9.0585E-01 

8.8729E-01 

83%8E-01 

030 

9.4707E-01 

9.2458E-01 

9.0223E-01 

8.8035E-01 

83934E-01 

8.2849E-01 

035 

9.2682E-01 

9.0132E-01 

8.7601E-01 

83133E-01 

8.2777E-01 

7.9376E-01 

0.60 

9.0403E-01 

8.7528E-01 

8.4663E-01 

8.1867E-01 

7.9216E-01 

7.5477E-01 

0.65 

8.7989E-01 

8.4721E-01 

8.1428E-01 

7.8199E-01 

73156E-01 

7.1005E-01 

0.70 

83613E-01 

8.1811E-01 

7.7903E-01 

7.4036E-01 

7.0414E-01 

6.5693E-01 

0.75 

8.3508E-01 

7.8909E-01 

7.4056E-01 

6.9193E-01 

6.4669E-01 

5.9097E-01 

0.80 

8.1971E-01 

7.6130E-01 

6.9783E-01 

6.3342E-01 

5.7398E-01 

5.0524E-01 

0.85 

8.1365E-01 

73568E-01 

6.4869E-01 

53953E-01 

4.7794E-01 

3.8934E-01 

0.90 

8.2120E-01 

7.1274E-01 

5.8933E-01 

4.6207E-01 

3.4665E-01 

2.2831E-01 

0.95 

8.4728E-01 

6.9223E-01 

5.1361E-01 

33904E-01 

1.6311E-01 

1.1938E-03 

1.00 
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Table  9:  Error  coefficients  Ff^„(a,  P)  -  pp  /brn  =  O(zeroth),  Ksecond), 

4(fou.rth),  6(sixth)  and  8(eighth)  for  a  coaxially  displace  coil  pair 


a  =  1.01,  P  =  05,  p,  =  035 


rla^ 

zeroth 

second 

fourth 

sixth 

eighth 

0.05 

1.4646E-03 

-1.1488E-06 

-7.3922E-09 

3.1827E-11 

-6.2614E-14 

0.10 

1.4646E-03 

-45953E-06 

-1.1828E-07 

2.0369E-09 

-1.6029E-11 

0.15 

1.4646E-03 

-1.0339E-05 

-S.9877E-07 

23202E-08 

-4.1081E-10 

0.20 

1.4646E-03 

-1.8381E-05 

-1.8924E-06 

13036E-07 

-4.1035E-09 

0.25 

1.4646E-03 

-2.8721E-05 

-4.6202E-06 

4.9730E-07 

-2.4459E-08 

030 

1.4646E-03 

-4.1358E-05 

-9.5803E-06 

1.4849E-06 

-1.0517E-07 

0.35 

1.4646E-03 

-5.6292E-05 

-1.7749E-05 

3.7444E-06 

-3.6096E-O7 

0.40 

1.4646E-03 

-7.3525E-05 

-3.0279E-05 

83433E-06 

-1.0505E-06 

0.45 

1.4646E-03 

-9.3055E-05 

-4.8501E-05 

1.6914E-05 

-2.6953E-06 

030 

1.4646E-03 

-1.1488E-04 

-7.3922E-05 

3.1827E-05 

-6.2614E-06 

035 

1.4646E-03 

-1.3901E-04 

-1.0ffi3E-04 

5.6384E-05 

-1.3422E-05 

0.60 

1.4646E-03 

-1.6543E-04 

-1.5329E-04 

93305E-05 

-16923E-05 

0.65 

1.4646E-03 

-1.9415E-04 

-2.1113E-04 

13362E-04 

-5.1076E-05 

0.70 

1.4646E-03 

-2.2S17E-04 

-18398E-04 

2.3964E-04 

-9.2405E-05 

0.75 

1.4646E-03 

-2.5848E-04 

-3.7423E-04 

3.6253E-04 

-1.6047E-04 

0.80 

1.4646E-03 

-Z9410E-04 

-4.8446E-04 

5.3397E-04 

-2.6893E-04 

0.85 

1.4646E-03 

-3.3201E-04 

-6.1741E-04 

7.6823E-04 

-4.3678E-04 

0.90 

1.4646E-03 

-3.7222E-04 

-7.7601E-04 

1.0825E-03 

-6.9001E-04 

0.95 

1.4646E-03 

-4.1472E-04 

-9.6337E-04 

1.4973E-03 

-1.0634E-03 

1.00 

1.4646E-03 

-45953E-04 

-1.1828E-03 

2.0369E4)3 

-1.6029E-03 

a  =  1.01,  p 

=  0.55,  pj  =  0.40 

r/ai 

zeroth 

second 

fourth 

sixth 

eighth 

0.05 

1.3858E-03 

-4.0654E-07 

-8.7867E-09 

2.8198E-11 

-3.5557E-14 

0.10 

1.3858E-03 

-1.6258E-06 

-1.4059E-07 

1.8047E-09 

-9.1026E-12 

0.15 

1.3858E-03 

-3.6580E-06 

-7.1172E-07 

2.0556E-08 

-2.3329E-10 

0.20 

1.3858E-03 

-6.5031E-06 

-12494E-06 

1.1550E-07 

-2.3303E-09 

0.25 

1.3858E-03 

-1.0161E-05 

-5.4917E-06 

4.4059E-07 

-1.3889E-08 

0.30 

1.3858E-03 

-1.4632E-05 

-1.1388E-05 

1.3156E-06 

-5.9722E-08 

0.35 

1.3858E-03 

-1.9916E-05 

-2.1097E-05 

3.3175E-06 

-2.0498E-07 

0.40 

1.3858E-03 

-2.6013E-05 

-3.5990E-05 

7.3919E-06 

-5.9665E-07 

0.45 

1.3858E-03 

-3.2922E-05 

-5.7650E-05 

1.4986E-05 

-1.5306E-06 

030 

13858E-03 

-4.0645E-05 

-8.7867E-05 

2.8198E-05 

-3.5557E-06 

035 

1.3858E-03 

-4.9180E-05 

-1.2865E-04 

4.9954E-05 

-7.6219E-06 

0.60 

1.3858E-03 

-5.8528E-05 

-1.8220E-04 

8.4199E-05 

-1.5289E-05 

0.65 

1.3858E-03 

-6.8689E-05 

-2.5096E-04 

1.3611E-04 

-2.9005E-05 

0.70 

13858E-03 

-7.9663E-05 

-3.3755E-04 

2.1232E-04 

-5.2475E-05 

0.75 

13858E-03 

-9.1450E-05 

-4.4483E-04 

33119E-04 

-9.1128E-05 

0.80 

13858E-03 

-1.0405E-04 

-5.7585E-04 

4.7308E-04 

-1.5272E-04 

0.85 

13858E-03 

-1.1746E-04 

-7.3388E-04 

6.8063E-04 

-14804E-04 

0.90 

13858E-03 

-1.3169E-04 

-9.2240E-04 

93908E-04 

-3.9184E-04 

0.95 

13858E-03 

-1.4673E-04 

-1.1451E-03 

1.3266E-03 

-6.0388E-04 

1.00 

1.3858E-03 

-1.6258E-04 

-1.4059E-03 

1.8047E-03 

-9.1026E-04 
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3.7  Axial  Field  Flux  in  the  Midplane  Between  a  Pair  of 
Coils 


From  Figures  24  and  25,  it  is  dear  that  if  the  axial  field  HJir/a-i,  n/2)  was 
integrated,  then  there  would  be  a  partial  caitcellation  of  fields  greater  than  Hg 
with  those  less  than  Hg.  Ultimately,  in  the  course  of  calibrating  an  induction 
coil  placed  between  a  pair  of  coils,  the  required  experimental  quantity  is  the 
flux  /  B  dA.  Referring  to  Figure  24  for  example,  one  could  choose  smaller 
values  of  <  0.3)  and  maintain  an  averaged  uniformity  over  a  laiger  area 
by  increasing  r/a-i.  However,  this  leads  to  undesirable  large  positive  and 
negative  variations  in  HJx/a-i,  jt/2)/Hg  from  unity.  Furthermore,  beyond  a 
certain  r/a-i  value,  the  field  deviation  drops  off  very  quickly  and  good 
uniformity  cannot  be  achieved. 

Because  of  cylindrical  symmetry,  it  is  possible  to  obtain  a  simple  expression 
for  the  flux  of  the  axial  field  in  the  mid-plane  of  the  coil  gap, 

=  }ig2jc  /  Hj(^,  n/2)4  d^.  Here  the  zonal  harmonic  expansion  form  is 
convenient  because  the  field  geometry  factors.  Table  2,  are  independent  of  4. 
The  flux  <]>  normalized  to  that  of  a  uniform  field  ((l>g)  obtains  from  equations  33 
and  38, 

[fg  -.|[foE2(a,P) 

+  ^[FgE,(a,P)  -F^E,(a^p^)]  -^[FgE,(a,p)  -F^E,(a^pp]^^°^ 


_  1 


For  a  =  1.01,  p  =  0.5  and  P^  =  0.35  (Fig.  24),  H^(0.65, 7t/2)  differs  by  3%  from  Hg, 
but  the  flux  within  the  central  mid-plane  region  bounded  by  a  circle  of  radius 
0.65  a-i  differs  by  only  0.5%  from  the  flux  of  a  homogeneous  field.  Typical 
values  of  {(r/flj),  «I>/Oj  evaluated  for  the  above  coil  geometry  from  eqn.  40  are 
as  follows:  (03),  1.0062;  (0.4),  1.0095;  (0.5),  1.0114;  (0.5726)^  1.0103;  (0.6), 
1.0091;  (0.65),  1.0050  and  (0.7),  0.9979.  For  a  =  1.01,  p  =  0.55  and  p^  =  0.4,  coil 
geometry  (Fig.  25),  typicjd  values  of  {(r/a,),  ®/«l>g)  are:  (0.2),  1.0010;  (0.3), 
1.0015;  (0.4),  1.0010;  (0.45),  0.9998;  (0.5),  0.9977  and  (0.6),  0.9887. 


’  (H.(r/«„  n/2)/Ho)  -  100 


4.  Helmholtz  Coil  Apparatus 


4.1  Introduction 

Based  on  the  computations  given  in  section  3,  a  Helmholtz  coil  pair  was 
constructed  to  provide  good  field  uniformity  over  radial  distances  of 
approximately  one  metre.  Since  the  lengths  of  the  coils  which  require 
c^brating  will  be  short,  the  field  homogeneity  Hjizlai,  OI/Hq  above  and  below 
the  central  plane  of  the  coil  pair  is  not  critical.  Furthermore,  in  order  to 
calibrate  smaller  coils  with  radii  significantly  less  than  one  metre,  it  was 
decided  to  have  a  second  coil  gap  setting  such  that  the  field  homogeneity  could 
be  improved. 


4.1.1  Physical  Description  of  Coil  Pair 

A  Helmholtz  coil  pair  was  constructed  in  order  to  achieve  the  field 
homogeneity  predicted  in  Figures  24  and  25,  specifically  a  =  1.01,  ^  =  0.50, 

=  0.35  and  a  =  1.01,  P  =  0J5,  P^  =  0.40  where  a,  =  99  cm  and  consequently 
02  =  99.99  a  100  cm.  This  leads  to  an  overall  coil  length  (including  gap)  and 
gap  spacing  of  99  cm,  693  cm  and  108.9  cm,  79.2  cm  respectively  for  the  above 
mentioned  P,  P^  values.  The  cross-section  in  each  coil  is  a  rectangle  1  cm  by 
14.85  cm. 

The  2  m  diameter  coil  pair  needed  to  be  placed  in  a  stable  frame  to  ensure  the 
correct  coil  spacing  while  allowing  easy  access  for  coils  which  require 
calibrating.  The  framework  consisted  of  four  stout  vertical  wooden  posts 
which  supported  two  square  wooden  platforms.  These  in  turn  secured  the 
coils,  maintained  the  correct  coil  gap  and  held  the  vertical  posts  in  a  rigid 
framework.  The  entire  structure  was  made  from  timber,  glue  and  nylon  bolts. 
Within  the  coil  gap,  two  horizontal  stout  timber  beams  were  cross-braced 
diagonally  to  the  vertical  posts.  This  structure  which  was  vertically  adjustable 
formed  the  support  for  the  coil  which  was  being  calibrated. 

The  two  square  platforms  which  support  the  coils  contain  wooden  dowel  pegs 
at  20®  spadngs  in  the  form  of  a  circle.  The  pegs  have  two  functions;  firstly  to 
secure  the  coils  in  the  correct  position  such  that  the  coil  pair  is  concentrically 
aligned  along  the  vertical  axis.  Secondly,  the  two  square  platforms  when 
bolted  together  make  up  the  demountable  former  upon  which  both  coils  were 
wound.  Here  the  pegs  hold  the  cheeks  apart  at  the  correct  spacing  and  also 
interlock  at  10®  spacings  to  provide  the  pseudo  circular  support  for  a  plastic 
strip  which  forms  the  mandrel.  This  technique  enables  both  coils  to  be  wound 
on  the  same  coil  former  thus  helping  to  preserve  the  same  geometry. 

A  horizontal  coil  pair  configuration,  rather  than  a  vertical  one,  enables  a  more 
even  weight  distribution  of  the  coils^  and  also  allows  the  coil  which  requires 
calibrating  to  be  more  easily  aligned.  In  order  to  minimize  weight,  aluminium 
wire  was  chosen  instead  of  copper  because  its  density  is  less  than  one  third  of 


^  The  colls  would  not  be  able  to  support  their  own  weight  in  a  vertical  configuration  and  would 
thus  require  further  structural  support. 


copper.  An  alternative  would  be  to  use  plastic  insulated  wire.  This  has  the 
advantage  of  requiring  fewer  turns  of  wire  for  the  same  a,  P  parameters,  the 
only  difference  being  the  larger  space  packing  factor  associated  with  the  plastic 
coated  wire  compared  to  enamelled  wire.  It  is  worth  noting  that  for  plastic 
insulated  wire  the  diameter  tolerances  of  the  conducting  filaments  can  be  as  low 
as  7%  according  to  Cable  Makers  of  Australia.  This  error  might  be  expected  to 
lead  to  an  overall  error  in  because  of  the  dependence  of  and  X  on  the 
conductor  cross  sectional  area.  The  homogeireity  is  unaffected.  However 
determines  the  overall  current  density  NI/l2b(a2  -  Ui)],  and  the  error  reduces  to 
that  of  2b{a2  -  a^).  Furthermore,  the  small  deviation  from  uniform  current 
density  distribution  caused  by  the  deviation  in  wire  diameter  would  be  small 
over  the  large  dimensions  of  the  coils. 

The  coils  were  woimd  with  2.00  nun  diameter  grade  2  polyester  imide 
enamelled  alununium  wire  having  the  following  properties:  = 

2.824  X  10^  Q  m  at  20‘’C,  -  2  7  x  10^  kg  m-^,  R  =  9  x  10^  ft  m-\  4.7  turns  per 

cm  and  20  tun\s/cm^.  This  implies  X  =  0.63,  =  16.8  kg,  R  >=>  17.80  ft  and 

N  B  300  turns.  These  figures  refer  to  each  individual  coil.^ 

4.1.2  Coil  Inductance  and  Mutual  Inductance 

In  order  to  minimize  the  inter-winding  capacitance,  inductance  and  coil 
winding  effort,  thick  wire  was  used  iiwtead  of  thin  wire.  The  inductance  (H) 
for  an  individual  coil  was  computed  from  Wheelers'  formula  for  short  coils  [8] 


L 


2.916  X  lO'^N^rlogjo 


4.97 

(flj  -  a,  +  2b') 


(41) 


which  gives  0.48  H  for  the  above  coils  (N  =  333)  and  compares  favourably  with 
the  measured  inductance  of  0.47  H.  In  eqn.  41,  r  is  the  mean  coil  radius, 

(flj  +  «i)/2.  The  mutual  inductance  M  was  computed  from  expressions  with 
tabulated  coefficients  representative  of  coil  geometry  [9]  and  gave  M  =  0.0696  H 
for  coils  of  rectangular  cross-sections.  (The  same  result  obtains  with  the  two 
coils  each  represented  by  a  single  circular  filament  multiplied  by  the  square  of 
the  number  of  turns.)  By  measuring  the  inductance  of  the  coil  pair,  both  series 
aiding  and  series  opposing,  the  mutual  inductance  was  determined  to  be 
0.068  ±  0.016  H. 


®  The  two  colls  were  individually  wound  with  five  layers  and  a  total  number  of  turns  N  =  333. 
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4.2  Theoretical  Fields 


4.2.1  Field  at  Centre  of  a  Non-zero  Cross-section  Coil 

As  an  approximation,  a  coil  can  be  replaced  by  a  single  filamentary  loop  where 
the  loop  radius  a  is  (02  +  «i)/2.  From  eqn.  2  at  N  =  333, 

Bj(0, 0)  =  [2.103  X  10^  J]  T.  For  a  coil  of  non-zero  cross-section  with  the  above 
parameters,  eqn.  7  leads  to  8^(0, 0)  =  [2.118  x  Iff^  J]  T  where  a  =  (100/99)  = 
1.0101.  The  single  filamentaiy  loop  coil  underestimates  the  central  field  by  less 
than  1%. 


4.2.2  Field  at  Centre  of  Coaxially  Displaced  Coils 

From  eqn.  33,  the  field  at  the  centre  of  the  gap  for  a  pair  of  coils  is  given  by 
B^(0,  0)  =  [3.284  X  IfT*  i]  T  for  a  =  1.01,  p  =  0.50,  =  0.35  and  B,(0,  0)  = 

[3.107  X  10^  1]  T  for  a  =  1.01,  P  =  0.55,  p^  =  0.40.  The  central  field  of  the 
Helmholtz  coil  is  55%  or  47%  higher  th^  the  single  coil  depending  on  the  gap. 


4.3  Determination  of  the  Calibrating  Field 

4.3.1  Search  Coil 

The  field  at  the  centre  of  the  coil  pair  can  be  predicted  from  the  expressions  for 
Bj(0,  0)  in  section  4.2.2  from  the  current  in  the  coils  which  can  be  monitored 
from  the  voltage  aaoss  a  1  Q  resistor  in  series  with  the  coils.  The  resonant 
frequency  of  the  coil  pair  and  the  individual  coils  was  measured  at  about 

2.1  kHz.  Thus  the  predicted  field  would  not  be  expected  to  be  accurate  above 
frequencies  where  the  inter-winding  capacitance  shunts  the  current  across  the 
coil.  In  this  case  the  voltage  drop  across  the  1  £2  sense  resistor  does  not  give  a 
true  indication  of  the  current  in  the  windings  and  hence  of  the  magnetic  field. 
Furthermore,  introducing  a  coil  at  the  centre  of  the  gap  can  further  load  the 
Helmholtz  coils  at  and  above  the  resonant  frequency  of  the  centre  coil.  Even  if 
the  centre  coil  is  open  circuit,  it  acts  as  a  shorted  turn  at  the  coil  centre  above 
resonance. 

A  convenient  method  of  determining  the  calibration  field  is  to  use  a  single 
layer  search  coil  wound  on  a  non  metallic  tube.  The  voltage  induced  by  a 
sinusoidal  magnetic  field  at  frequency  /,  orthogonal  to  the  plane  of  the  search 
coil  windings  is 


V  =  INv^r^fB 


(42) 


where  r  is  the  radius  of  the  winding  in  metres,  B  is  in  tesla  and  V  is  in  volts. 
One  hundred  turns  of  38  S.W.G.  enamelled  copper  wire  was  woimd  on  a  PVC 
tube  having  an  outside  diameter  of  155.5  mm.  The  search  coil  calibration 
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factor  is  11.93  V  Hz'^  or  equivalently 


B 


8.381  X  10 


-2  V 

T 


(43) 


The  resonant  frequency,  75  MHz,  of  the  search  coil  lies  well  above  the 
Helmholtz  coil  resonance.  (A  less  sensitive  search  coil  wound  with  a  hundred 
turns  (38  S.W.G.)  on  25.4  nun  diameter  rod,  having  a  calibration  factor  of 
03184  V  T’  Hz-^  was  also  used.) 

43.2  Comparison  of  Search  Coil  Measurements  and  Theoretical 
Expressions 

In  the  following  discussion,  refers  to  the  field  determined  from 
measurements  tal<en  with  the  more  sensitive  search  coil  (eqn.  43)  and  Bf^  refers 
to  measurements  taken  with  the  smaller  less  sensitive  sean^  coil.  In  addition, 
Bjq  refers  to  the  field  at  the  coil  pair  centre  obtained  from  the  theoretical 
expression,  8,(0,  0)  =  3.284  x  10^  I T,  by  measuring  the  voltage  drop  across  the 
one  ohm  resistor  in  series  with  the  coil  pair  windings.  All  voltages  were 
measured  with  a  spectrum  analyser.^ 

Below  about  5(K)  Hz,  determinations  of  the  central  field,  with  search  coils  and 
from  the  theoretical  expressions  gave  agreenient  within  4%  or  better.  The 
discrepancy  at  higher  frequencies  is  illustrated  in  Figure  26  which  plots  Bf^/Bi^ 
as  a  function  of  frequency.  Similar  results  obtain  for  B/VBjq.  Qearly,  above 
500  to  600  Hz,  the  determination  of  the  central  field  from  the  measurement  of 
the  Helmholtz  coil  current  is  inaccurate.  This  frequency  range  is  further 
reduced  if  the  central  coU  which  requires  calibration  has  a  resonance  frequency 
below  the  Helmholtz  coil  pair  resonance  frequency. 


4.33  Comparisons  with  Fluxgate  Magnetometers 

The  field  from  the  Helmholtz  coil  pair  was  also  measured  with  a  Domain 
magnetic  underwater  sensor  (MUWSl  triaxial  fluxgate  system)  and  a  less 
sensitive  fluxgate  probe  (Institut  Dr  Forster  Magnotoscop  1.068).  The  field 
measured  by  these  two  sensors  are  designed  B,  and  Bp  respectively.  The  only 
significant  restriction  with  these  two  instruments  is  that  they  are  limited  to  low 
frequency  operation,  especially  the  MUWSl.  Measurements  were  taken  at 
1  Hz  and  05  Hz  over  a  range  of  coil  currents  with  care  being  taken  not  to 
saturate  the  MUWSl  sensor.  The  estimated  errors  are  B/VBjq  :  ±  3%, 

Bp/BjQ  :  ±  3%,  and  B,/Bjq  :  ±  2%.  Good  agreement  to  within  experimental 
error  was  found  between  Bjq  and  B,  and  between  Bjq  and  Bp.  Reasonably 
good  agreement  was  also  found  for  Bj”/B,q  (1.043  ±  0.020). 


^  In  this  way  one  instniment  allowed  measurement  of  voltages  over  the  extremely  low  frequency 
range  and  furthermore,  the  waveform  could  be  monitored  to  allow  for  distortion  or  non-linearities. 
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Frequency,  Hz 

Figure  26:  Ratio  of  magnetic  field  measured  with  search  cod  (Bf^)  to  that  measured 
from  current  in  Helmholtz  cod  (Bjq)  as  a  function  of  frequency.  (At  high  frequencies, 
B]q  is  not  derioed  from  the  current  in  the  cod). 

The  above  measurements  were  testricted  to  /  ^  1.0  Hz.  At  higher  frequencies 
up  to  100  Hz,  the  search  coil  field  measurements  were  checked  with  a  known 
calibrated  field  generated  by  the  vertical  component  of  the  three  axis  cube 
arrangement  of  the  magnetic  facility  at  the  Shark  Point  degaussing  range, 
Sydney.®  A  series  of  magnetic  field  measurements  with  both  search  coils 
reveal  an  agreement  generally  within  ±  2%  of  the  computed  field  at  the  centre 
of  the  cube. 


’  Here,  agreement  to  within  1%  was  found  with  die  MUWSl  at  05  Hz. 
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4.4  Field  Homogeneity 


The  field  homogeneity  of  the  coil  pair  is  readily  determined  by  taking  field 
measurements  using  a  search  coil  at  known  distances  along  one  of  tlw  arms  of 
the  cross-beams,  see  4.1.1.  The  search  coil  windings  are  located  in  the  central 
plane  between  the  Helmholtz  coil  pair.  Measurements  were  also  taken  for 
displacements  above  the  centre  of  the  coil  pair.  Thus  the  theoretical  curves  of 
HJir/a-i,  ic/2)/Ho  and  HJiz/a-i,  0)/Hq  in  Figure  24  can  be  verified  experimentally. 

The  results  of  the  field  uniformity  tests  are  given  in  Table  10.  They  show 
excellent  agreement  with  the  theoretical  curves. 


Table  10:  Magnetic  field  homogeneity  assuming  a  =  1.01,  P  =  0.5,  P^  =  0.35, 
Figure  24 


H,  (t  n/2)/Ho 

Experimental  Theoretical  (Table  3.2) 

1.0 

1.0 

0.0 

1.0 

1.002 

0.1 

1.0 

1.006 

0.2 

1.012 

1.011 

03 

1.012 

1.015 

0.4 

1.012 

1.012 

03 

0.977 

0.992 

0.6 

0.923 

0.936 

0.7 

0.813 

0.812 

0.8 

0.603 

0369 

0.9 

H,  (Y,  0)/Ho 

1 

Experimental 

Theoretical  (Table  3.2) 

1.0 

1.0 

0.0 

0.994 

0.997 

0.1 

0.977 

0.986 

0.2 

0.966 

0.966 

0.3 

5.  Conclusion 


The  design  for  air-core  uniform  current  density  coils  is  discussed  in  terms  of 
coil  geometry.  This  enables  the  field  at  the  centre  of  the  coil  (Hq)  to  be 
computed  or  determined  graphically  from  Figures  3  and  4. 

The  magnetic  field  within  the  region  of  the  coil  centre  is  evaluated  using 
Garrett's  [3]  zonal  hormonic  expansion  technique.  This  enables  the  radial  and 
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axial  fields  to  be  expressed  as  a  siim  of  terms  involving  geometry  factors  and 
normalized  distances.  The  direct  dependence  of  the  field  expressions  on  Hg 
allows  the  fidd  homogeneity  to  be  conveniently  expressed  in  terms  of  geometry 
factors  alone.  Suitable  convergence  of  the  series  expansion  must  be  achieved. 
Table  3  gives  the  individual  contributions  of  the  second-,  fourth-,  sixth-  and 
eighth-order  terms  over  a  wide  range  of  coil  geometries.  Several  important 
exan^les  are  given  in  Figures  6  to  8. 

Hart's  [2]  method  for  computing  the  radial  and  axial  fields  is  also  discussed. 
This  involves  a  direct  evaluation  of  the  fidd  by  numerical  integration  of  a 
volume  dement  within  the  coil.  A  conq?arison  of  this  method  is  made  with 
the  zonal  harmoiuc  expansion  approach.  The  former  technique  is  a  useful 
con^utatioiud  method  for  an  in^vidual  coil  geometry  but,  unlike  the  expansion 
method,  is  not  suited  for  modelling  studies. 

The  superposition  principle  aiul  symmetry  argtunents  are  next  employed  to 
compute  the  axial  fields  along  the  axis  of  two  coaxially  displaced  coils.  A 
comparison  of  the  field  homogeneity  generated  by  this  system  with  that  of  a 
single  coil  reveals  the  improvement  offered  by  a  two-coil  system. 

The  effectiveness  of  the  zonal  harmonic  expansion  method  becomes  more 
apparent  when  it  is  applied  to  the  multiple  coil  system  to  provide  expressions 
for  the  fields  and  their  homogeneity  within  the  central  region  of  the  gap  formed 
by  the  pair  of  coils.  The  correct  gap  width  for  maximum  hontogeneity  or  for 
perfectly  compensated  coils  can  be  obtained  dther  graphically  (Figs  19  to  22)  or 
numerically.  The  graphical  approach  gives  a  good  approximate  result  and 
clearly  shews  the  correct  gap  setting  as  a  function  of  coil  geometry.  A 
partic^ar  result  that  emerges  from  Figures  19  and  21  is  that  by  choosing  a 
value  of  P  such  that  F  E-jst,  P)  is  a  minimiun,  then  for  small  a,  the  correct 
Helmholtz  coil  gap  equals  the  internal  radius  of  short  coils,  »>  p.  The 
deviations  associated  with  this  rule  of  thumb  for  larger  values  of  a  can  also  be 
gauged  from  Figures  19  to  22. 

The  numerical  solution  of  the  correct  gap  setting  for  fourth-order 
compensated  coils  allows  the  contributions  of  the  remaining  higher  order  terms 
to  be  automatically  determined.  This  allows  the  computation  and  tabulation  of 
error  limits  for  different  coil  geometries  (Table  6). 

Equations  32  and  33  allow  modelling  of  the  field  homogeneity  for  a  wide 
range  of  gap  settings.  Fourth-order  compensation  is  included  by  substituting 
the  requi^  P^  value,  see  Figures  19  to  22  and  section  3.4.  The  resultant  field 
homogeneities  HJir/a-i,  k/D/Hq  are  plotted  so  as  to  examine  the  trade-off 
obtained  by  achieving  a  laiger  area  of  field  uniformity  at  the  expeitse  of 
degrading  the  degree  of  field  uniformity.  Depending  upon  the  application,  this 
can  be  advantageous  since  the  total  flux  of  the  axial  field  in  the  midplane  can 
approach  the  value  that  would  exist  for  a  perfectly  uniform  field  because  of  the 
partial  cancellation  of  the  positive  and  negative  areas  enclosed  by  the 
Hj(r/aj,  Jc/2)/Ho  versus  r/ui  cxirves.  This  advantage  does  not  apply  to  the 
monotonically  decreasing  function  0)/Hq. 

In  section  4,  a  physical  description  and  construction  details  are  given  for  a  coil 
pair  based  on  the  theoretical  expressions  developed  in  section  3.  The  field 
predictions  were  checked  with  the  aid  of  two  search  coils  and  two  fiuxgate 
magnetometers.  Agreement  between  theory  and  experiment  lies  within  4%. 
Furthermore,  the  field  homogeneity  was  found  to  be  in  excellent  agreement 
with  theoretical  expectations. 
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Appendix  A 

Numerical  Computations 


Fortran  source  code  (HELMHZ.for)  was  developed  for  evaluating  the  axial  field 
deviation  from  the  central  field  for  a  Helmholtz  coil.  The  percentage  deviation 
from  the  central  field  is  read  in  as  the  variable  ERROR.  Ibe  variable  DELTA 
determines  the  tolerance  in  the  Newton-Raphson  solution  for  determining 
Sample  computations  cue  shown  in  Table  6. 

Fortran  source  code  (HSPUT.for)  was  also  written  to  compute  the  deviation 
of  the  axial  magnetic  field  from  the  central  field  at  points  located  along  the  axis, 
Hj(z/ai,  0)/Ho  and  radially  away  from  the  gap  centre, 
refers  to  the  axial  field  at  the  centre  of  the  gap.  Computations  are  carried  out 
over  a  range  of  specified  values  for  a  given  a,  P  value;  a,  P  and  are 
defined  m  Figure  A.I.  (Uidike  Figure  14,  2fi  is  equivalent  to  the  overM  coil 
length  including  the  gap,  whereas  previously  in  Figure  14,  2^  referred  to  the 
length  of  each  individual  coil). 


Figure  A.1;  Coaxial  coil  pair. 


The  following  variables  are  required  from  a  data  file. 


AC,  BC  the  a,  P  values 

GMIN  :  the  initial  value  of  p/a,  (z/a,  or  r/a,) 

NGAMMA  :  the  number  of  increments  in  p /a, 

DELG  the  p/a,  increment 

N6GAP  the  number  of  values 

BETA(I)  the  values 


A  subroutine  COEFS  computes  the  required  error  coefficients  f(a,  P),  f  £2(0,  P) 
....fEg(a,  P)  in  double  precision. 

These  programmes  may  be  requested  from  the  author  if  desired. 
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